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-.^ , Recently, gauge field theory approaches were extensively used in order to discuss the physical conse- 

■ quences of spin-orbit interactions in condensed matter physics. An SU(2)xU(l) gauge theory is very 

naturally borne out and provides an illustrative example of a classical Yang-Mills field theory at work. 
This approach may serve as an exemplification of non-Abehan field theories for students in general 
physics curriculum. It allows to introduce discussions on fundamental ideas like Noether currents, gauge 
symmetry principle, gauge sjonmetry breaking and non linear Yang Mills equations in very concrete 
physical situations that makes it accessible to a broad audience. 

Resume. Les theories de jauge ont recemment ete utilisees de maniere systematique pour discuter 
les consequences physiques de I'interaction spin-orbite en physique de la matiere condensee. Une theorie 
SU(2)x U(l) apparait naturellement et fournit un exemple d'illustration du fonctionnement des theories 
classiques de Yang et Mills. Cette approche pent servir d'example de theorie de jauge non abelienne 
pour des etudiants dans un cours de physique generale. Cela permet d'introduire et de discuter les idees 
fondamentales comme les courants conserves de Noether, le principe de symetrie de jauge, la brisure de la 
symetrie de jauge et les equations non lineaires de Yang et Mills dans des situations physiques concretes, 
les rendant ainsi accessible a une vaste communaute. 

ip{r, t) wave function or complex scalar field 

i/)(r, t) Pauli 2— component spinor or Pauli spinor field 

t) Dirac 4— component spinor 

if Lagrangian density 

H Schrodinger Hamiltonian 

H Pauli Hamiltonian 

®t, ® U(l) covariant derivatives 

<p, A U(l) scalar and vector gauge fields 

£, B electric and magnetic fields 

p(r, t) charge density 

j(r, t) or J(r, t) charge current density 

e, g electric and isospin charge 

T° SU(2) generators 

a SU(2) contraction T° a" 

St, @ SU(2) covariant derivatives 

s" a— component of the spin i operator 

<T° a- component Pauli matrix 

= |(T°;f °, W = ^cr^y^l" SU(2) scalar and vector gauge fields 

E = ifT^E", B = i<T°B" SU(2) field strengths 

^mat f)' Pmat f*^' ■^"(r, t) a— Component of the matter spin or isospin density and current density 

^rad Prad </*0^j a— Component of the radiation spin or isospin density and current density 

str, t) vector spin density 



I. INTRODUCTION 



Yang-Mills theories are usually introduced to students in the context of quantum field theory, but are never encountered 
in their classical form at the undergraduate level fT] . Regarding the fundamental role of gauge symmetry in the construc- 
tion of physical theories, this is desirable in order to make the underl5ang concepts accessible to physics students early in 
the curriculum. In the abstract of a paper by R. Mills [2], the content of gauge theories is summarized as follows: ". . . the 
gradual emergence of symmetry as a driving force in the shaping of physical theory; the elevation of Noether's theorem, relating 
symmetries to conservation laws, to a fundamental principle of nature; and the force idea ("the gauge principle") that the sym- 
metries of nature, like the interactions themselves, should be local in character " This obviously brings the concept of gauge 
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invariance to the level of a keystone principle among physical theories, rather than a specific property of electromagnetic 
theory. 

The non relativistic limit of Dirac equation, namely the Pauli equation and its corrections to Schrodinger theory, i.e. the 
spin-orbit (SO) interaction and the Zeeman interaction, offers a natural frame for an SU(2) gauge field theory. This does 
not only give an elegant formulation for the coupling of spin degrees of freedom with external electromagnetic fields, but 
also the Yang-Mills approach provides natural answers to fundamental issues such as conservation laws which may offer 
illuminating insights in the condensed matter physics context |Q] • Here we show that it may be a useful and very concrete 
way to introduce Yang-Mills theory, in the condensed matter physics curriculum. 

A reformulation of the spin-orbit coupling Hamiltonian in terms of non-Abelian gauge fields [0-01 was explicitly given 
in ref. where the SO interaction is presented as a SU(2)x U(l) gauge theory. A key property for that purpose is to 

deal with an SO interaction which is linear in the momentum, H50 ~ p • (<? x £) in such a way that it can be absorbed in the 
kinetic energy through a minimal coupling procedure, p — > p — x £, up to unessential numerical factors. This property 
is fulfilled by the Pauli Hamiltonian, but also by the Rashba SO interaction [fl4tl and by the two-dimensional reduction of 
the Dresselhaus SO interaction [fisll . A potential realization of such situation in the general case is thus the 2— dimensional 
electron gas (2DEG) as encountered in semiconductor hetero-junctions [fT6ll and 2DEG quantum rings |tT^. 

Such gauge point of view, in more general terms, has been known for some time[|18l-[20[l. This formulation is very 
revealing, since the consistent gauge structure of the theory becomes obvious and the physics of spin currents, persistent 
currents and color diamagnetism[21] can be understood in a manner analogous to the well known U(l) gauge theories. 

The plan of the paper is as follows: we start in section [ll] with the case of electromagnetic theory which will then serve 
as a basis to the non-Abelian Yang-Mills theory presented in section [nil Section [W| discusses, in this context, the case of 
spin-orbit and Zeeman interactions and application to Einstein-de Haas effect is discussed in section |Vl We have decided 
to avoid most of the time the use of tensor analysis and instead embrace vector algebra, in order to obtain equations of 
motion in a form very similar to Maxwell equations in matter, as they are known in standard textbooks [i2a1 . We feel 
that this will help the part of the community who is less familiar with field theory, appealing to the broader intuition from 
electromagnetism. 



II. A WELL-KNOWN GAUGE FIELD THEORY: THE EXAMPLE OF ELECTROMAGNETISM 

In this section, we collect the standard results of classical gauge field theory applied to electromagnetism in order to 
later follow a parallel approach in the case of non Abelian and SO interaction. 

A. Conserved currents 

A fundamental question in physics poses the problem of the electric charge. Why is charge conserved? Historically we 
have made the empirical observation that in all fundamental processes, e.g. collisions between elementary particles, the 
electric charge is conserved. The sum of the charges of the outgoing particles is always equal to the sum of the charges of 
the incoming particles. This property then translates into a hydrodynamic conservation equation d^p{Jr, t) + V -j^f, t) = 0, 
where p(r, t) is the electric charge density calculated on a coarsed grained scale when a continuous density make sense. 
Integration over a mesoscopic (or macroscopic) scale, then leads to = J d^r d^p{r, t) = j{r, t)- dS = provided 
that the volume of integration is large enough and the current densities vanish at the boundaries. So conservation of 
charge follows from an experimental observation, but its origin is not something obvious. An answer to the question of 
the origin of such a conservation law (formulated here in a provocative manner) is that Hilbert spaces and the complex 
fields living in these spaces don't associate physical relevance to global phases of such fields. Here we refer to the fact 
that in Quantum Mechanics (QM) the results of measurements are given by the matrix elements {Lp\Obs.\ip) of Hermitian 
operators Obs. associated to physical observables, and such matrix elements are independent of the choice of reference for 
the phase of the complex vector c^. 

The principle of conservation of the electric charge, experimentally legitimate, is derived using the second approach 
from a more fundamental and more abstract principle of invariance of the theory with respect to global phase changes. 
The main advantage of the second answer is that it enables us to connect to other similar theories. More importantly, it 
leads to a consistent construction of the electromagnetic theory with no a priori knowledge other than that implied by the 
Lagrangian formalism which is the language adapted for such approaches. This is what we now show. 
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The Schrodinger equation follows from a least action principle [[2 
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where the free matter Lagrangian density ifg is an ordinary function ^ of the wave function ip{f,t) (a complex scalar 
field describing matter), its space and time derivatives, and the corresponding complex conjugates treated as independent 
variables. ££q = ^{^, V(/?, Vip*, d^ip, d^ip*) is the Jordan-Wigner Lagrangian density [l23l.[24fl . 
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and equation ([T]) applied to this Lagrangian recovers the Schrodinger equation. As discussed above. Quantum Mechanics 
is invariant under a global change of phase of the wave function 



c^(r, t) —> exp{iea/h)ip(r, t), a = const. 



(3) 



since all physical quantities, expressed in terms of matrix elements, are independent of the global phase choice a. The 
coefficient e/h in the phase is present for later convenience. It simply fixes the dimensions of the ea to those of an action. 
The phase transformation above is a continuous s)Tnmetry (a can take any real value) which leaves the Lagrangian density 
unchanged. Applying Noether theorem to an infinitesimal gauge transformation ip —> ip' [l + iea/fl)(p = ip + Stp [l23l.[25[l . 
it follows that we can define a conserved current obeying a continuity equation. Indeed, demanding the global gauge 
transformation to be a symmetry, i.e. to preserve the lagrangian density ifg unchanged leads to 



5ifo = ■■ 



5ifo d^o ^ 

-5y + . 5(p -\ =; — 5V(^ + ^ ) 



d(p 



dVip 



lea 



d, ip 



dip* dip 



ip +V- ip 



— - — V 

dVip 



-a{d,p + V-j) 



(4) 



where use has been made of the equations of motion ([T]l in the second line. This leads to the definition of a charge 
density p{r, t) = eif*{f, t)ip{j, t) and a charge current density )(r, t) = t)V(/j(r, t) — ip{r, t)Vip*{r, t)). This 

very interesting connection between "phase invariance" in QM and charge conservation was first made by F. London. For 
interested readers, historical aspects of gauge invariance are discussed at length in the book of O'Raifeartaigh [f2^ or in 
the reviews [f27-.2al. 



B. The gauge invariance principle illustrated 

Let us now quote Salam and Ward [29], (see also Novaes [|30|]1, who speculated on whether a local generalization of 
the gauge invariance would generate strong, weak and electromagnetic interactions and the associated new fundamental 
conservation laws: "Our basic postulate is that it should be possible to generate strong, weak, and electromagnetic interaction 
terms ... by making local gauge transformations on the kinetic-energy terms in the free Lagrangian for all particles." The 
message behind this sentence is that theories which have a larger symmetry class (local gauge symmetry) have richer 
dynamical structure (new interactions present), and the power of gauge field theory appears when one assumes that the 
phase invariance invoked above also holds in the case of local phase changes a(r, t), 

ipir, t) exp(iea(r, t)lh)ip{f, t) = Gip{r, t), (5) 

where G = G[j,t). This is obviously not a sjmimetry of the original Lagrangian density ^q, since Vt^(r, f) and dfip{f,t) do 
not transform like the wave function itself, and as a consequence, terms like ihip^d^ip or ~^^'^*'^'^ ^re not gauge invari- 
ant (the first of these terms for example transforms into ihip*diip — e[d^a)ip*ip). In order to repair this lack of symmetry of 
the original Lagrangian, and to be able to build a new one with the expected invariance property, one may define so-called 
covariant derivatives @ = V — ^A(r, t) and = 3^ + ^<^(r, t) in such a way that acting on ip{f, t), they transform like the 



4 



wave function, i.e. according to 2)if[r, t) — » exp(iea(r, t)/h)^(fi[r, t) and ®t(/'(r, t) — » exp(iea(r, t)/h)^i-(f{r, t). The two 
fields A(r, t) and (f){r, t) are introduced for the transformations above to be satisfied. Again, the factors e/h are present 
for convenience and the imaginary i is a choice which enables and A to be real fields (we require id^ and — iV to be 
Hermitian operators). The gauge transformation conditions imply that these fields obey the laws 

Air,^ GAlir,t')G-^ + —i\/G)G-^ =Alir,t) + Va(r,t') (6) 

le 

4>ir,t) - G0(r,OG-i--(a,G)G-i = 0(r,O-5,a(r,t), (7) 

le 

where one recognizes the gauge tranformations of the vector and scalar potentials of electromagnetism. 

Summarizing the discussion to this point, we made a local phase change to the wave function, but this does not leave 
the Lagrangian density invariant unless we compensate this change i.e. the terms brought about by derivating a can be 
cancelled by corresponding changes in the gauge fields. The terminology of covariant derivative refers to the fact that 
and are not invariant under a gauge transformation, but do transform in the same manner as the field if itself. 
The derivative operators and @ themselve transform like vectors under a gauge transformation, G&^G'^ and 

^ ^ G&G~^ . In the presence of an external gauge field, one way to make the Lagrangian density gauge invariant is the 
minimal coupling prescription, i.e. replace the ordinary derivatives by the covariant ones in the free matter field Lagrangian 
density 

= J^o(i^,^*,...) + ^i^,i^,^*,...,4>,A), (8) 

where now there appears a supplementary interaction term, depending on both the "matter field" and the gauge fields 
(f), A. This term appears because any term allowed by the gauge symmetry should appear In principle higher order terms 
can also appear but one takes the minimal lowest orders. The price to pay to extend the so-called U(l) gauge invariance 
to local phase transformations is that it automatically generates the interactions of the charged matter field with the 
electromagnetic field, which is a strong self consistent check for the theory! Yang summarized the gauge invariance 
approach to the study of fundamental interactions by the expression "Symmetry dictates interaction" [jsil] 

In terms of the gauge fields, the physical electric and magnetic fields E and B are defined as the commutators of 
covariant derivatives, £^ = ^[®t)®yt] and = — |^eijjc[®,-, ®j] (eijjc is the totally antisymmetric tensor which takes 
value +1 when ijk is 123 up to circular permutations, it is —1 when ijk is 213 up to circular permutations, and it is 
when two or three indices are identical). In 4— dimensional notation, one introduces the Faraday tensor F^,, = ^ [@^, @^.] 
(with @^ = 3^ + ^A^ and the usual metric signature, A^ = {(p /c, —A)) which transforms under gauge transformation as 
a vector, i.e. F^^ GF^^.G~^. The Faraday tensor (hence the electric and magnetic fields) is thus gauge invariant in 
the Abelian U(l) gauge theory. The contribution of the free gauge fields to the Lagrangian density then follows simply 
from the following requirements: we demand a i) Lorentz invariant, ii) quadratic in the gauge fields derivatives (i.e. a 
kinetic energy of a standard form) and iii) gauge invariant quantity. Prescription i) is fulfilled if all space-time indices 
are contracted while ii) is nicely obtained from the square of the Faraday tensor The quantity f^^f which transforms 
under gauge transformation as F^^f^^ ~* GF^yf^^ G~^ satisfies all three conditions, since G and G~^ simplify each other 
in the present Abelian theory for which all quantities commute. This free gauge field Lagrangian density is known as the 
Schwarzschild Lagrangian density [|32l - [34l] ). 

^m) = -i-/,rF'''' = ISoE'-i-B\ (9) 

The formulation of the principle of gauge invariance was stated in Yang and Mills work [Q], where they generalize 
"phase invariance" of ordinary Quantum Mechanics of charged fields. They always refer to the case of electromagnetic 
theory in order to illustrate their arguments: We define isotopic gauge as an arbitrary way of choosing the orientation of 
the isotopic spin axes at all space-time points, in analogy with the electromagnetic gauge which represents an arbitraiy way 
of choosing the complex phase factor of a charged field at all space-time points. We then propose that all physical processes 
(not involving the electromagnetic field) be invariant under the isotopic gauge transformation ip —>■ ip', ip' = S'^ip, where S 
represents a space-time dependent isotopic spin rotation. 

To preserve invariance one notices that in electrodynamics it is necessary to counteract the variation of a [the a of Yang and 
Mills is the same as our —ea in Eq. ([5])] with x, y, z and t by introducing the electromagnetic field A^ which changes under 
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a gauge transformation as 

In an entirely similar manner one introduces a B field in the case of the isotopic gauge transformation to counteract the 
dependence ofS on x, y, z and t. (. . .) The field equations satisfied by the twelve independent components of the B field, which 
we shall call the b field, and their interaction with any field having an isotopic spin are essentially fixed, in much the same 
way that the free electromagnetic field and its interactions with charged fields are essentially determined by the requirement of 
gauge invariance. 

It was also clearly stated by Utiyama [@] : The form of the interactions between some well known fields can be determined by 
postulating invariance under a certain group of transformations. For example, let us consider the electromagnetic interaction 
of a charged field Q(^x), Q*(x). The electromagnetic interaction appears in the Lagrangian through the expressions d^Q — 
ieA^Q, d^Ct + ieA^Ct (!)• The gauge invariance of this system is easily verified in virtue of the combinations of Q, Q*, and 
A^^ in (1), if this system is invariant under the phase transformation Q —> e'°Q, Q* Q*e~'° a = const (2) [because the 
combination — ieA^Q also transforms as Q]. Reversing the argument, the combination (1) can be uniquely introduced 
by the following line of reasoning. In the first place, let us suppose that the Lagrangian L{Cl,Q ^) is invariant under the 
constant phase transformation (2). Let us replace this phase transformation with the wider one (gauge transformation) 
having the phase factor a(x) instead of the constant a. In order to make the Lagrangian still invariant under this wider 
transformation it is necessary to introduce the electromagnetic field through the combination (1). This combination and the 
transformation character ofA^ under the gauge transformation can be uniquely determined from the gauge invariance postulate 
of the Lagrangian L(^Q, Q,^,A^). 



C. Conserved currents and equations of motion in the presence of gauge fields 

In the present context of an Abelian gauge theory, the free gauge field contribution ifu(i) does not play any role in the 
expression of the conserved currents. This is due to the fact that the bosons (photons here), which are responsible for the 
(electromagnetic) interaction, do not carry the corresponding (electric) charge. We shall see later that this is no longer true 
in the non Abelian case. Nevertheless, the conserved current gets an additional term due to the presence of the interaction 
of the matter field with the gauge potentials [Issll . Indeed, since any local phase transformation is now a S3mimetry of 
the Lagrangian if'^at. = + -^int.- this is also true for global phase transformations and equation Q now becomes 

5(^o+^mt.) = = -a(5,p + V-J) (10) 

where the charge density p(r, t) is unchanged, but the current density acquires an additional term often called diamagnetic 
current density, J(r, t) = ^iv*ir, t)^^pif, t) - t^(r, t\§if{r, t))*) = /(r, t) - (e^/m)A(r, t)if *{r, t)ip[i^, t). In terms 
of this conserved current density, the interaction term is often written to lowest powers in the gauge fields as if'int = 
J - A— pep + 0(^), where all fields depend in the general case on the variables (r, t). 
The total Langrangian density becomes 

-^tot. = -^mat. + -^Lf(l) 

+ieo(-V0-3,A)2-5i-(VxA)2 (11) 

The equations of motion follow from Euler-Lagrange equations (O. For the matter field equation of motion, we perform 
the variation w.nt. t^*(r, t). 
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ihdtLp = {-ihV -eAfip + e(pi^ + Vip. (13) 
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Those for the gauge fields also obey Euler-Lagrange equations, but follow from variations w.rt. the gauge fields, 
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and the very definition of the charge and current densities p 



■ g^"- , J = g J'- . The Maxwell equations follow 



V X (V X A) - eoMo^tC-Vc/) - S,A) = /XqJ, 



V-(-V</)-3,A) = -, 



(14) 
(15) 



(16) 
(17) 



the two equations without sources being automatically satisfied with the definition of the physical fields E and B in terms 
of the gauge fields. In the presence of other types of interactions of the matter field, there can occur new contributions to 
the charge and current densities in Maxwell equations as we will see later. 

The gauge field theory approach thus appears as a very powerful tool in order to define properly the notion of conserved 
currents, the interaction of a matter field with an external (gauge) field, the equations of motion of these gauge fields, as 
well as some other subtle notions such as topological invariants |^] . This last notion is very useful for example in the 
description of the Aharonov-Bohm effect in which a charged particle moves in a pure gauge field (i.e. a region of space 
where the physical field B vanishes but the potential vector does not). In such a case, the gauge invariant phase of the wave 
function ^^Adr is responsible for the shift, proportional to the magnetic flux enclosed by the path, of the interference 
pattern produced by any kind of two-slit experiment [IsT 



III. YANG-MILLS THEORY 



Our program is now to follow the same lines of derivation as previously outlined, but in the case of a SU(2) symmetry. 
This is done in many excellent textbooks [;38:-43:], so we will just motivate the theory and fix the notations that we will 
need in the application we have in mind to condensed matter physics. The U(l) gauge transformation introduced in the 
previous section was motivated by the fact that non relativistic (spinless) quantum mechanics deals with wave functions 
which are complex scalar fields and that physics is independent of the choice of global phase of this field. A complex scalar 
field can also be seen as two real scalar fields (the real and imaginary parts), and the gauge transformation thus appears 
as a rotation in the complex plane [l44t] . We can imagine that a more elaborate theory would require the introduction of 
more sophisticated mathematical objects, e.g. spinors ||45l] (or a set of complex scalar fields). In the simplest case, we 
will assume that the state of the particle is no longer described by a simple complex scalar field, but by a Pauli spinor 

consisting in two complex scalar fields usually written in the form of a "vector" i/; = and we will call SU(2) (special 

unitary group of 2 x 2 matrices with determinant 1) the corresponding gauge theory. This name comes from the fact that 
the objects acting on a two component spinor are two by two matrices which can all be written in terms of four basic 
ones, the unit matrix and the Pauli matrices. Now, we obviously have a set of 4 real scalar fields, but only 3 are in fact 
independent because the wave function must be normalized as J d'^r i^'^-i^ < +oo where the adjoint spinor is the "row 
vector" ip'' = (ifil (/^p. In the "internal" space called isospin space (name given by Eugene Wigner) where these three 
components live, a gauge transformation may be seen now as a rotation in a 3— dimensional space. Since the action of 
SU(2) operations rotate in three dimensional space and such rotations do not commute (order of rotations changes results) 
as in two dimensions, one says that the theory is non Abelian. This contrasts with the Abelian U(l) theory where complex 
scalars commute under multiphcation. Such a transformation, acting on the spinor, may be written as 

i/>(r, t) exp(iga/h)xp(r, t). (18) 

The constant g plays the role of a new charge and all operators acting on the spinors are defined as 2 x 2 matrices (as 
expected), denoted here with bold font to emphasize their matrix nature, e.g. a stands for the contraction over the 
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repeated index a r'^a'^ =T^a} + T^a^ + t^o^ where the r" are the generators of SU(2) (i.e. half the PauH matrices), 

The a^'s are the parameters of the transformation, they define the components of the rotation vector in the three- 
dimensional isospin space. 

The transformation ( [18] ) is a symmetry when the parameters a^'s are just constants, hence a conservation equation 
follows from Noether theorem, since, as quoted by R. Mills, "The substance of the theorem, for our purposes, is that for every 
symmetry of nature there is a corresponding conservation law and for every conservation law there is a symmetry" and its 
extension to local gauge symmetry is the strategy originally proposed by Yang and Mills [4^ : "The conservation of isotopic 
spin points to the existence of a fundamental invariance law similar to the conservation of electric charge. In the latter case, 
the electric charge serves as a source of electromagnetic field. An important concept in this case is gauge invariance which is 
closely connected with (1) the equation of motion of the electromagnetic field, (2) the existence of a current density, and (3) 
the possible interactions between a charged field and the electromagnetic field. We have tried to generalize this concept of gauge 
invariance to apply to isotopic spin conservation". Isotopic spin or (strong) isospin had been introduced earlier by Heisenberg 
in 1933 as a new conserved physical quantity in order to explain the close similarity of behaviour of protons and neutrons 
during strong interactions, e.g. the same cross section (apart from electromagnetic corrections) for interactions involving 
both t5^es of particles, like p + n~ n + or n + tz* ^ p + n°. The nucleon was then considered as an isospin doublet. 

The extension of the transformation (flST ) to local parameters a°(r, t), faces the problem that exactly like in the case of 
U(l) symmetry, the space and time derivatives of the spinor do not obey the same transformation law ( flSl) as the spinor 
itself does. As we mentioned before the simple gradient applied to the spinor is not gauge invariant, and thus is not a 
valid building block for the Lagrangian. The alternative is to come up with a recipe for a new derivative, the covariant 
derivative, which after having acted on a spinor, behaves also like a spinor (as one should expect, the derivative of a spinor 
is another spinor). 

A space covariant derivative has to be defined as 

# = VI - f\N, (20) 
or in component form ©j. = dj^l— ^r^W", and a covariant time derivative is similarly introduced by, 

&r = d,l + fx = 3,l + fT^X"- (21) 

1 is the 2x2 identity matrix. The non-Abelian gauge potentials and carry an isospin index, and they obey specific 
transformation laws through local gauge transformations. 



[/(I) Global gauge transformations Si7(2) 

ihip*ip — ^(V(/7)*(V(/?) free particle Lagrangian density ££q ihili 'ip — ^(Vli/^J'XVli/)) 

ihdtifi = —^V^ip free particle equation of motion ihd^lip = — ^V^lip 

if — > exp{iea/h)ip gauge transf. ip —> exp{^ga.)il) 

dfP+V-j = continuity equation S^p" + V • j° = 

p = e(f*ip charge density p° = gi/i'Vi/i 

j = ^(^^p*'^^p-^p^^p*) conserved current /" = ^(i/'V''(V'i/>) - (V'i/')V°t/j) 

Table 1: Comparison between U(l) and SU(2) global gauge invariance. 
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Once the non-Abelian gauge potentials have been defined, we can form a field tensor and a gauge covariant Lagrangian 
density in order to allow for a simple derivation of the conserved current density following from Noether theorem, and 
essentially follow similar lines of derivation than in the previous section. In order not to repeat all calculations, we present 
the analogies between U(l) and SU(2) theories in two tables. The properties associated to global gauge invariance are 
listed in the first table, while the consequences of local gauge invariance follow in the next table. 

Let us briefly comment upon the results listed in these tables. In Table 1 there is a very clear one to one correspondence 
between the quantities defined in the two gauge theories. When one invokes global gauge invariance in the SU(2) case, 
isospin is inert fixed in a certain direction and everything operates as if charge is the only physical quantity. This is 
consistent with the interpretation of the Schrodinger equation not as an equation of motion for spinless particles, but 
for "two-component Pauli particles" in a spin eigenstate [l47ll . There are nevertheless important differences between the 
Abelian and the non-Abelian theories. First, the fact that the gauge potentials carry internal indices (this is the origin of 
the non-Abelian character), implies that the same indices are carried by the conserved charges, as well as by all other fields 
such as E° and B°, the non-Abelian analogues to the electric and magnetic fields defined in Table 2. Another essential 
difference is due to the fact that the gauge fields also carry the "charge" of the interaction. This has two main consequences 
(see Table 2) : the conserved charge density and current density have, together with the matter part, a purely radiative (i.e. 
involving only the gauge fields) contribution and the equations of motion for the gauge fields are non-linear 



1/(1) Local gauge ^^(2) 

transformations 



if ^ exp{iea{r,t)/h)ip = Gip local gauge transf. t/) ^ exp(^gct(r, t))!/" = Gi/' 

f = V-fA 1 ( @ = V1-^W 

" COV. derivatives ^ " 



St = 3, + '-§^ I [ = 3,1+ 'fx 

A^A+Va \ gauge transf. | W ^ W Va + f [a, W] 

>P^<p-dtaj of potentials y ^ ^ ^ - g^ct+ f[a,x] 

^^ = ^^^^] 1 field strengths h^ = UfM 

p = eip*ip cov. charge density Pmat. +Prad. = S'^^ ^''^' 

J = ^(V*(Vy) - y(V(^.)*) cov. conserved J" + T = ' -r^CVV-) " (V^-)' t°^) 

mtW = ^{-ihV -tAf^ + ti,^ matterfield ihdA.->h = g'Hf-^ + gvtb 

^ equation of motion ' ^ 2m^ 

V X B - c-^d,E = iiJ ] gauge fields | V x B° - c-^S.E" = ^(J" + /°) 

V.£=f^ J equations of motion | ^ . = i-(p-^^^+ p-^J 



Table 2: Comparison between U(13 and SU(2) local gauge invariance. 



Another subtle difference between the Abelian and the non Abelian case resides in which quantities are gauge invariant, 
which are gauge covariant, and which are simply gauge dependent. One is used to the case of the electromagnetic 
interaction for which the charge density and the total current density are gauge invariant (even though in the latter 
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case the paramagnetic and diamagnetic contributions are separetaly not gauge invariant, their sum is). The electric and 
magnetic fields (or the Faraday tensor) are similarly gauge invariant and finally only the wave ftmction is gauge covariant 
ip' = G(f. The covariant derivative operator introduced to satisfy the same covariance law (©^v')' = Si'^'f' = G'2iy_^p 

then tranforms as a vector = G@^G~^ and this implies that the gauge potential, although a 4— vector, has a particular 
non-covariant gauge transformation = GA^G~^ — j^{d^G}G~^ . In the non Abelian situation, the same consideration 
applies for the spinor \p' = Gip and the covariant derivative operator, = G@^G~^, and this requires the gauge potential 
to obey a similar non-covariant gauge transformation W^^ = GW^G~^ — -^(5^G)G~^. But now the so-called curvature 
tensor F^,, transforms like a vector, F^^, = GF^G~^, i.e. is subject to a rotation in the isospin space, hence the SU(2) 
magnetic and electric fields are not gauge invariant and are then unphysical in a gauge symmetric theory (for a discussion 
on the meaning of "physical meaning", see e.g. Ref. [l4^ ). The non Abelian charge density and current density themselves 
are not gauge invariant, and under a gauge transformation, their matter and radiation content changes. Only the total 
non Abelian charge Q" = J d^r p'^{r, t) is gauge invariant (and constant in time) provided that the volume of integration 
extends over large enough distances for the current to vanish. 



IV YANG-MILLS THEORY FOR SPIN-ORBIT INTERACTION 



A. Preamble - the Zeeman interaction 



Electrons are particles which carry not only an electric charge, but also a spin |. In non-relativistic Quantum Mechanics, 
they obey the Pauli equation which is a generalization of the Schrodinger equation with an essential innovation: the 
existence of the spin, i.e. a new degree of freedom, coupled to space-time degrees of freedom in the Hamiltonian through 
interactions such as the spin-orbit interaction and the Zeeman interaction. In the presence of such interactions, spin is 
known not to be conserved, but spin carries angular momentum and the total angular momentum is conserved. This is at 
the origin e.g. of the Einstein - de Haas experiment that will be discussed later. 

Although these statements are well understood physically, a clear formulation is appealing following the lines of deriva- 
tion proposed in the quotation of Yang and Mills |46|] in section|lIll Here the starting point is given by the known interaction 
terms present in the Pauli equation, from which a consistent gauge theory is built. Such a formulation leads, as before, to a 
conservation law. Then, once the gauge transformations are extended to local changes, it leads to the equations of motion 
for the gauge fields. 

In this introductory section, we will focus attention on the consequences of the Zeeman interaction between the magnetic 
moment associated to the spin of the electron (with Lande factor ge=2 and s = and an external magnetic field 

B. The state of the particle is described by a Pauli spinor ip = { ), and the Hamiltonian reads as H = r'-Cp — eA)^l -I- 

VI — j^s- B. Following standard textbooks [|4^.l5ol] one can build a continuity equation where the charge density and the 
charge current density are defined according to 

p(r,t) = ea/)'i (r, t)i/.(r, t), (22) 

J(r, t) = ^ixP\n t)[VV;(r, t)] - [^xp\r, t)]V(r, t)) 
2m 

e^^ e ^ 

Aip^ + —V X [xP^'sip). (23) 

m m 

As we had already anticipated, when a new interaction (here the Zeeman interaction) appears in the problem, the con- 
served current gets modified. This is a property that we have already seen at play when the minimal coupling to the 
external electromagnetic field was added to the free particle, where it led to the appearance of a diamagnetic contribution 
to the Noether conserved current density. Here, the rotor of the magnetization associated to the electron density produces 
an additional charge current density. Note that this last term may be forgotten during the standard derivation using the 
wave equation and its complex conjugate, since the divergence of a rotor vanishes, and some care must be taken to estab- 
lish the full current density [jsil - lssl ]. Nevertheless, this addititonal term, sometimes called "spin term" [52], is compulsory 
in order to obtain a conservation equation. It is thus also instructive to contemplate the spin density (e.g. in the direction 
of the external field B = Be^), which may be naturally defined by 5^(r, t) = i/^ ' s^i/). The z— component commutes with 
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the Hamiltonian. This is a conserved quantity, and the corresponding continuity equation takes the form 5^5^ + V ■ = 
with J" = ^[ip 's^^xp — ['V\p'^)s^\p) the spin current density. This is a tensorial quantity, since it has a vector character 
(the current propagates in space) and depends also on the spin orientation. The conservation equation also follows from 
the use of Ehrenfest theorem, ^ = [ih)~^ip'' [s^, H]ip + ip'^d^s^'ip = 0. The other components of the spin however do not 
commute with the Hamiltonian and in this case the continuity equation gets additional terms. Altogether, one can write 
3^5° + V • (s^v) = ^(sx B)a, where the RHS represents the torque exerted on the spin. 

The spin is just a contribution to the total angular momentum and this latter quantity is of course conserved. Our 
discussion thus suggests that there should be a way to write a conservation equation where spin components would appear 
explicitly, together with other sources of angular momentum. The "spin density" components s°(r, t) = ip^'s^xp are, in a 
sense equivalent to the charge density, except that they carry the spin index a. Deriving such a conservation equation is 
the purpose of the next sections. 



B. Pauli equation and spin-orbit interaction 

We are now in a position to discuss completely the role of the spin-dependent relativistic corrections to the Schrodinger 
equation. The Pauli equation follows from the non relativistic limit of the Dirac equation [ImII . 

ihd,^ = [cd ■ip-eK) + l3mc^ + V] *, (24) 

where a and 13 are the 4x4 Dirac matrices and * is a 4— component Dirac spinor. We let * = ^ ^ j g-imc^t/n remove the 

rest energy of the electron, xj) and £, both have two components. In the non relativistic limit, ^ is a very small component 
f =i — J — C(T • (p — e^xl), and xh satisfies 

md^lip=H'4). (25) 

Here, a is the vector of Pauli matrices a = 0-^6^^ + a-^Cy + cr'e^ = cr'^e^. Forgetting about the relativistic correction to the 
kinetic energy, the Pauli Hamiltonian is given by 

H = — (p - eA)^! +V1 
2m ^ 



Non Relativistic KE 

eh eh ^ ^ ie}i^ _» _ eh: 

- ^^-^ -ir^^-^Ex{p-eA))-—^a-TIxE + —^V^. (26) 
2m 4m^c^ 8m^c^ 8m^c^ 



Zeeman term SO interaction Darwin term 

In the following, we will neglect the Darwin interaction which, like the relativistic correction to the kinetic energy, is 
spin-independent. The "substrate" potential V will be denoted as ecp, defining the scalar gauge potential. Let us note that 
the SO interaction is written in its more general form here, including the term proportional to the rotor of the electric 
field [i5J0- 



C. Gauge theory of the spin-orbit interaction 

The fact that the SO interaction is linear in p enables one to gather it together with the Non Relativistic KE in a single 
squared term, ^ [(p — eA)l — gW]^ up to a correction quadratic in the spin operator. Since icr" is the generator of SU(2) 

sjmimetry (denoted as r° in tables 1 and 2), we set W = icr^W", with W" three ordinary vectors. Expanding the cross 
products in terms of antisymmetric tensors, the SO interaction in Eq. ([26] ) may now be written as 

Hso = -^^i^"eajkEjPk + '-^(^"e^jkdjE^'] (27) 
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and identified to the relevant contribution of the cross term in the Non Relativistic KE, 

Hso = -^a''(p,W^+W^p,) = ^a''(2W^3,+(a,W^)), (28) 

provided that the new "SO" charge, here denoted as g, is identified to g = ft. It is thus natural to define a new gauge vector 
potential, depending on both space-time variables and spin degrees of freedom, gW = ^^<? x E [I'ssI] or in components 

Iga^W" ^ igcT^^Nle, = ^^e,,,a%e,. (29) 

The same identification is obtained if we include the electromagnetic vector potential. We see that the gauge vector 
components carry two indices. The index k here refers to a purely space index while the index a is a spin index which 
also has a spatial meaning due to the fact that spin and space variables are linked in Quantum Mechanics (the index a 
would be a purely internal isospin variable in the context of Yang-Mills theories of fundamental interactions). In order to 
avoid any confusion between the two indices, we will always denote as upperscript the spin indices (the position of the 
indices makes no difference, since we use a Euclidean metric for this non relativistic theory). Repeated indices are summed 
up, independently of their up or down position. The SO interaction is now essentially included in the KE through 

Non Relativistic KE -I- SO Interaction = 

1 , ^ eft e^ft^ 
[pl-eAl \^^E\^. (30) 

Let us point out the important fact that the gauge vector potential has an inherent matrix structure, hence the underlying 
gauge theory is a non-Abelian one, called SU(2) due to the presence of the Pauli matrices, exactly like what was briefly 
presented in section [ini The addititional term in equation dsUj) has the structure of a gauge symmetry breaking term [|57l] 
i.e. it changes when we perform a gauge transformation. This, we now know how to do for this theory, and its effect will 
appear in a later discussion. 

In the case of the SO-l-Zeeman interactions of equation (|26| ). one has to build a full U(l)xSU(2) theory. This means 
that we contemplate the space of complex scalar and their products and the space of 2 x 2 matrices and their product. 
Both these spaces being independent. To take into account simultaneously the fact that the electric and magnetic fields are 
minimally coupled to the charge and to the spin degrees of freedom, we define the covariant derivatives 

® = VI - ^^1 - iW, 

(31) 

where A is the electromagnetic vector potential and W° = ^^e^ji^Ep while 4> is the usual electromagnetic scalar potential, 
and 2° = — — B°. Note that we have used the simplification which occurs due to the fact that the non-Abelian charge is 
g=ft. 

Under an infinitesimal local gauge transformation G = exp(iicr''a°) 1 -I- i^cr'^a", the non Abelian gauge fields should 
transform according to 

X^X' = GXG-' - -(5,G)G-i = la^ix" - e,,,a'x' - d.a^ (32) 
W^W = GWG-^-F-(VG)G-i = icr''(W''-e„f„a''W'^ + Va''). (33) 



The alert student should note that in this theory of the spin-orbit interaction there is no gauge freedom associated to 
SU(2) transformations (actually a small subset survives fiol]) since the non-Abelian gauge fields are defined in terms of 
the physical Maxwell fields [|55f]. We will come back later to this point. Another peculiarity is due to the appearance 
in equation (|30]) of the so-called gauge symmetry breaking term (GSB) (the last term in the equation) that interferes 
with gauge symmetry. The full Lagrangian hence becomes ^^^^ = ifg + -^int. + -^gsb + ^u[i) + ^suii)' where ifmat. = 
ifo + -S?int. + -^GSB (-^mat. coutaius all tcrms involving the matter field xp}. We also have to add the free field contribution 
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^sv{2) (called Yang-Mills Lagrangian density), which only depends on the non-Abehan fields. This term is required by 
gauge symmetry and Lorentz invariance (for the theory to be automatically extended to a relativistic context). 

^ V ' " V ' 

+ - c^B.Bj + A(E^E^ - c^B^B^), (34) 

^ V ' ^ V ' 

with 

P = -3,W«-Vz'' + e„b,x''W^ (35) 
B° = V X W" + ie^jb.W'' X (36) 

following from the definitions E;^ = ^[5^1 + ix,3{\. — iW^] and B;^ = —Y^eij]^[^^l — iW;,3jl — iW^] given in Table 2. 
The dimensions of the fields are consistent with the coefficients appearing in the covariant derivatives, [j"] = T~^, 
[W^] = [E^] = £~^T~^ and [B^] = £~^, where T and £ are time and space dimensions. The fact that the pure non- 
Abelian field contribution ^su(,2) is gauge invariant is not a priori obvious due to the gauge co-variance of the physical fields 
in the non-Abelian situation E;^ GE^^G"^ = E^ -I- i|a°[cr'',E;^] and the similar transformation law for B;^. Nevertheless, 
-^s[/(2) being proportional to TrF^^.F'^'' TrGF^,,F'^^'G~^, it is obviously gauge invariant thanks to the invariance properties 
of the trace under cyclic reordering of the matrices. There remains the question of the physical origin of this term in the 
present situation. As argued by Tokatly [|2lll . it has the correct form to account for the first SO corrections to the energy. 
In any case it is demanded by S3mimetry so it should be out there! 



D. Conserved currents 



Let us stress on the fact that, contrary to the results presented in Table 2 in the case of a SU(2) gauge symmetric 
theory (Lagrangian invariant under SU(2) gauge transformations), the so called dia-color contribution ——W^ip'ip to the 
current density is exactly compensated in the present situation by the GSB contribution. The term dia-color was coined 
in this context by Tokatly and refers to the analogue of the diamagnetic current density, but for the non Abelian SO 
interaction II21II . The matter current = finally simplifies from the expression in Table 2 to 

k 

-ih\ , 

K = ^ ' i^'^^i^^^ - (5,1/-)' i^"!/^), (37) 
while the radiation current = is the same as indicated in the table, 

k 

X = -'eabcix% + |c2e,,.,Wf Bp. (38) 
The GSB term has no effect on the matter and radiation spin densities which remain those of the table. 



8^, 



mat. 



mat. g 

and 



: xp''lh(T''%p (39) 



SLf(2) _ ft ...f, 



-^e,,,\NlEl. (40) 



The spin-dependent part in the interaction term if-jK is thus exactly J^W^ "■^mat ^° (with no quadratic contribution this 
time). A very simple, but important consequence is that there is no spin current contribution proportional to the non 
Abelian gauge vector, i.e. no current transverse to the ordinary electric field, and the spin Hall conductivity vanishes 
automatically in the GSB scenario [fioll . 
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E. Equations of motion 



In order to obtain the equations of motion for the gauge fields, we also have to take into account the last peculiarity we 
will consider of the present problem, i.e. the fact that the U(l) and the SU(2) gauge fields are not independent. There 
are two approaches we can follow to contemplate this dependence within the variational formulation: one either includes 
constraints using the Lagrange multipliers [fssi] method, which allows treating the different gauge fields as independent. 



dtd^r 



(41) 



or one has to include in the equations of motion an account of the second order derivatives in the U(l) gauge fields, e.g. 



5S 
5^ 



34) 



-3, 



tot. 



+d. 



g^tot. 

+ 23 A 



+ 3,3, 



3id,3,<p) ''3{d,dj4>) 



0. 



(42) 



Let us illustrate these two approaches in the very simple case of electrostatics where the Lagrangian density reduces to 
^ = |eo|£|^ — p0. Using the first route, the condition E = —Vcp is implemented in the action through the constraint 

S = J d^r [if — A ■ (£ + V(/))]. Varying the action w.rt. (f) (£ being considered as independent of 0) leads to |^ + ^ ■ 

A • (£ + V<^) = — p + V • A = 0, while variation w.r.t. E produces the equation EqE — A = 0. The combination of the 



two equations recovers the usual Maxwell equation. The second approach is more straightforward in this simple case and 



3^ 



—p — BqV cj) = 0, i.e. the Poisson equation. 



leads to If -V g^^^^ 

The first route is by far simpler in the SO case. For the usual electric and magnetic fields, we obtain the following 
equations of motion (Maxwell equations) : 



5S„ 



5S„ 



54> 



54) 



+ ■ 



5S 



L/(l) 



54> 



5if 



30 



5 

= ~Pmat. + ^O^kEk + -^^3^^aki^% 

= 0. 



(43) 



5Ai 

3^„ 



5A, 



5A, 5A, 



dtd^r 



3 if, 



1/(1) 



3Ai 



d£^uc,) e 3 e 3 



'3{3kAd '3{d,Ad m '5(3,A,)^ 



2mc2 '5(3tA,) 



= 0. 



(44) 



In the first of these expressions, we can define a polarization in terms of the Lagrange multipliers, 

e 



Cnt; A? 



"/c o 9 ^ aki 



(45) 



Note that the origin of the dielectric polarization has to be found in the fact that a moving magnetic moment /x creates 
a dipolar moment ~ v x /3 [,56] . This dielectric polarization is associated to bound charges and there appears then a 
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contribution —d^Pf^ to the total charge density. This polarization also contributes the second equation of motion through a 
term added to the ordinary current density, d^P^, as well as an additional term ei^^d^nf, with 

m'' = --/i^ (46) 
m 

describing "Amperian" currents (regular charge currents). 

The electric charge and current densities are modified by the SO+Zeeman terms and the usual Maxwell equations 
become: 

V-E=-(p-V-P), (47) 

V xB-eo/io3t£ = jL(o['/ + 5tP + V X m]. (48) 
The variation with respect to the non-Abelian gauge fields leads to 



5x" ~ Sx" 5x" Sx" 



dtdh^\x'' + -B''), 



5^" dx^ ^' 
0. (49) 



5W? 5Wf 5Wf 5Wf 



V-E'' = -(5l,+5^^^,. + ,.''), (51) 



= J- - ie^bc [I'EJ + c^e,,,W}B^J - ftce,,-,5,B^ + fS.E;' - , 

= 0. (50) 

We recover the radiation spin density 5°^^ = -2eQ(,£W^E^ and the radiation spin current density = -^eabclx^^t + 
ic^eij^W|B^] in terms of which the corresponding Yang-Mills-Maxwell equations become 

c 

V xB° -^d.E" = ^iJ" + T -X''). (52) 
nc 

The non Abelian equations without sources take the form (obtained directly from the definitions (|35T | and (|36]) ) 

V X P 5,6° = e„b,(V X (x^'W^) + iS^CW'' x W^)), 

From the Yang-Mills-Maxwell equations, one can form a continuity equation describing the conservation of the total 
angular momentum density. Taking the divergence of Eq. ( |52T l. and using Eq. JSTI) . one gets 

5.(Ct. +4d.) + V • (i° + /") = -3,^" + V • r . 

The Lagrange multipliers formally appear in the continuity equation, but physically, the l.h.s. contains all sources of angular 
momentum incorporated in the problem [IssI], i.e. from free electrons in the conduction band encoded in the spinor xp 
(the matter spin density), and from angular momentum transfered to the lattice via SO interaction (the radiation angular 
momentum density). So the conservation of angular momentum in the system reads as expected 

5.(Ct.+<ad.) + V-(i" + /'') = 0. (53) 

We note that in the case of a static and homogeneous problem, the Lagrange multipliers just take the simple form = 
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V EINSTEIN - DE HAAS EXPERIMENT REVISITED 

The Einstein-de Haas experiment is one of the famous experiments of the beginning of the XXth century based on 
gyromagnetic phenomena [5S-.60l1 . In this experiment, a cylinder made of a non magnetized ferromagnetic material is 
suspended to a torsion wire and an external magnetic field is applied along the cylinder's axis. The cylinder acquires a 
magnetization, because spins orient to the field lowering the energy of the system and at the same time, the cylinder starts 
a rotating motion. This effect became an efficient experimental method for the measurement of the gyromagnetic ratio of 
various materials and, for example, proved that magnetism of iron is essentially due to the spin degrees of freedom. 

This experiment is well known, and well understood from very fundamental principles, since it relies on the conservation 
of angular momentum. When the material gets magnetized, the individual magnetic moments of the electrons point in 
a common direction, hence a total angular momentum appears which has to be compensated by an opposite angular 
momentum carried by the whole sample in order to conserve the initially vanishing value of the total angular momentum 
of the sample. Although this mechanism is very clear and unquestionable, there is no simple microscopic explanation of 
how this conservation of angular momentum is at play at the atomic level. See however[6I] for a beautiful description of 
the dynamics of the effect in terms of elastic twist on the lattice induced by the rotation of localized spins (spin-rotation 
coupling) and the coupling to the anisotropy field. The latter anisotropy is in fact a consequence of the existence, at a 
deeper level, of spin-orbit interactions. The purpose of this section is to provide such a microscopic description. Our 
approach is based on the gauge field description of spin-orbit interactions. The interest of the gauge field theory here is to 
provide a natural way for the introduction of conserved currents carr5ang angular momentum. 

An external magnetic field B = Be^ is applied to the cylinder of ferromagnetic material initially non magnetized, sus- 
pended along its z— axis to a torsion wire (hence the appearence of a non-Abelian scalar field j'^ in the Maxwell-Yang-Mills 
approach). The material gets magnetized. The spins acquire a common orientation and the corresponding magnetiza- 
tion is called the matter spin polarization. The sample is subject to a torsion around the vertical axis, and the radiation 
counterpart is the angular momentum acquired by the lattice. We can understand the physical situation as the appearence 
of a matter spin polarization s^^^ since the individual localized spins get polarized along the direction of the external 

magnetic field. As a consequence, V ■ 7^ and the Yang-Mills electric field aquires a non vanishing component E^. The 
localized spins interact with the lattice ions through an SO interaction (responsible for example for crystalline anisotropy), 
hence some components of W" do not vanish, creating, together with /'^ a radiation contribution to the spin polarization 

Let us analyse the equations in order to see whether an equation of rotation of the lattice follows. Here, the lattice enters 
the problem via the SO interaction, i.e. via the electric field (the "true" internal electric field) which is "rigidly" associated 
to the lattice ions. Since s^^^ = — (ft/c)e3j,cW^E^, and E^ = —d^W^^ + e^^^x'^'^l assume here that = , determined 
by the applied magnetic field, is uniform) we obtain 

4d. = X + + (£3)^)5)3. (54) 

As we had noticed, the appearence of the radiation contribution to the polarization is connected to the SO interaction 
(hence the electric field). This equation self-consistently shows that the true electric field produced by the ions located 
at the lattice sites is subject to a motion of rotation around the axis 63 along which the external magnetic field was 
applied. This rotation of the electric field can only be caused by an equivalent rotation of the whole lattice, thus explaining 
the observed Einstein-de Haas effect. If we now take into account the magnetization acquired by the cylinder and the 
exchange interaction between the conducting electrons and the localized moments (J^^j coupling constant), B has to be 
changed into B^ff = B-^J^^\XqM. 

VI. CONCLUSION 

In this paper, we have shown that spin-orbit interaction is a suitable arena to introduce non-Abelian gauge field theory 
(namely here SU(2) gauge theory) and that many consequences, relevant in a condensed matter physics context can be 
drawn from the classical form of the theory. More specifically, we have discussed the conserved current associated to this 
theory, the equations of motion of the gauge field, and the role of the continuity equation in the discussion of the Einstein 
- de Haas experiment. The approach follows from a very general Yang-Mills-Maxwell theory which we believe is likely 
suitable for micromagnetic simulations of spin current and spin configurations in various materials such as conducting 
ferromagnets or 2DEGs. 
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